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A 

The  dynamical  behavior  of  laser-stimulated  surface  processes 
(LSSP)  is  studied  by  the  generalized  Langevin  equation  via  the 
memory  effects  of  the  damping  kernel,  the  dephasing  kernel  and 
mode-mode  interactions.  The  temperature-dependent  averaged 
energy  absorption  rate  (which  characterizes  the  line  shape)  is 
calculated  by  solving  the  Laplace-Fourier  transform  of  the 
velocity  correlation  function.  The  features  of  the  response 
function  and  its  overall  line  broadening  are  discussed  in  terms 
of  Markovian  processes.  The  nature  of  LSSP  and  laser-selective 
effects  are  discussed  in  terms  of  the  multiphonon  coupling  strength 
and  the  related  "internal  resonant"  condition.  It  is  shown  that 
laser-selective  bond  breaking  is  possible  for  a  slow  intramolecular 
vibrational  relaxation  rate  which  is  governed  by  a  high-order 
multiphonon  process  or  a  far  off  '‘internal  resonance.* 

A; 


I .  Introduction 

Two  general  problems  in  laser-stimulated  surface  processes 
(LSSP)  are:  (i)  the  dynamical  behavior  of  the  excited  adspecies/ 
surface  system  and  (ii)  the  selective  and  nonselective  nature  of 


LSSP . ^  ®  The  latter  problem  was  recently  studied  by  a  quantum 

model  where  the  energy  populations  of  a  multilevel  system  were 

numerically  analyzed  for  different  sets  of  the  pumping  rate,  the 

relaxation  rate  and  the  damping  factor. The  energy  absorption 

profiles  were  also  studied  by  the  classical  Langevin  theory,  where 

the  many-body  effects  of  the  bath  modes  are  effectively  replaced 

by  the  bath-induced  damping  factor  and  frequency  shift  in  the  Wigner- 
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Weisskopf-type  approximation.  In  the  present  paper,  we  shall  fur¬ 
ther  study  these  two  major  problems  of  LSSP  via  the  generalized 

Langevin  equation  (GLE)  and  the  multiphonon  coupling  strength  with 
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the  related  "internal  resonant"  condition.  ' 

The  relaxation  dynamics  of  LSSP  described  by  GLE  with  memory 
effects  is  presented  in  Section  II,  and  the  temperature-dependent 
averaged  energy  absorption  rate  for  the  model  systems  of  (i)  a 
single  active  mode  and  (ii)  two  coupled  active  modes  is  calculated 
in  Section  HI.  In  Section  IV,  we  analyze  the  nature  of  laser-selective 
effects  via  multiphonon  processes  and  derive  (by  a  classical  model) 
the  required  condition  for  a  rapid  intramolecular  relaxation.  The 
averaged  energy  absorption  rate  of  a  system  with  many  active  modes 
is  derived  (by  a  quantum  model)  in  the  Appendix. 


II .  Generalized  Langevin  Equation  and  Relaxation  Dynamics 

We  consider  a  heterogeneous  model  system  consisting  of 
species  (atoms  or  molecules)  chemisorbed  on  a  solid  surface  and 
subject  to  low-power  infrared  laser  radiation.  The  Lagrangian  of 
the  system  may  be  written  as 

j  £*#•$«?)  ♦ 

<■“■1 

where  (i=l , 2 , • • • ,N)  are  the  normal  coordinates  of  the  system 
(with  normal  mode  frequencies  .  f^ (t)  is  the  transformed  gener¬ 
alized  force  given  by 


£  t Qt  ~  Vw. (QiAr- X) 


(2) 


where  is  the  effective  electric  field  of  the  laser  radiation 

acting  on  the  i-th  normal  mode  (with  classical  charge  q^) . 

Vanh  ^l'^2' *  *  * the  anharmoni-c  mode-mode  coupling  potential 
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given  by 

V  ■  =  ex,  -  V  1 1 1  +  •  •  • .  *> 

L  1  •  j  L.  *  LL.  i  ■ 

where  X,.  and  X.  ^ 


V 


<P  . 

V 

are  the  coupling  constants  characterizing  the 


single-phonon  and  multiphonon  processes,  respectively.  The  equa¬ 
tions  of  motion, 

ft 

0; 


2  0; 


C4‘ 


are  a  set  of  coupled  nonlinear  equations.  As  pointed  out  in  a 

7 

previous  paper,  the  numerical  solution  of  the  above  equations 

involves  long  computation  times  due  to  the  high  frequency  of  the 

14  -1 

laser  field  (tu  -10  sec  )  ,  where  the  time  scale  of  the  energy 


absorption  profiles  is  much  longer  than  the  oscillation  cycle  of 
the  field  for  the  low-power  excitation  processes.  A  new  set  of 
transformed  equations  in  the  rotating  frame  was  developed  to  over- 


3 


7 

come  such  numerical  difficulty.  However,  for  a  three-dimensional 
solid  system,  the  low-frequency  bulk-phonon  modes  are  condensed 
and  give  us  essentially  an  infinite  number  of  normal  modes,  where 
any  numerical  method  will  be  unfavorable. 

For  tractable  results,  we  shall  model  the  system  in  the  fre¬ 
quency  domain  as  follows:  one  (or  more  than  one)  of  the  vibrational 
inodes  of  the  total  system,  which  strongly  interacts  with  the  laser 
field,  may  be  singled  out  as  the  pumped  active  (A)  mode,  and  the 
remaining  modes  are  divided  into  two  subsystems:  the  surface- 
phonon  (B)  modes  (with  high  frequencies)  and  the  bulk-phonon  (C) 
modes  (with  condensed  low  frequencies) .  The  condensed  C  modes  are 
considered  to  provide  a  heat  bath  for  the  (A+B)  modes.  By  employ¬ 
ing  the  Wigner-Weisskopf-type  approximation,  the  many -body  effects 
of  the  condensed  C  modes  may  be  replaced  by  introducing  a  damping 
factor  and  a  frequency-shift  in  the  A  and  B  modes.  Instead  of  the 


coupled  Equations  (4) ,  we  thus  may  describe  laser-stimulated  surface 


processes  by  the  generalized  Langevin  equation  (GLE)  as  follows: 

jit)  +  ^  QSt'l  +  or.  |  jS.  (i-i)  Q-tt')  H  '  +  rrf„  Q;  j  f/.  j-/'  %-A  J 1 

o  0 

.  (5' 


0 

0^  is  the  j-th  normal  coordinate  of  the  (A+B)  subsystem;  8j(h)  and 
Mj (t)  are  the  damping  kernel  and  the  dephasing  kernel,  respectively, 
describing  the  interaction  dynamics  among  the  (A+B)  modes  and  the 


bath  (C)  modes;  Nj,  j  (t)  is  the  interaction  function  between  the  i-th 
and  j-th  modes;  f j (t)  and  f j (t)  are  the  driving  force  of  the  laser 
field  and  the  bath-induced  random  force  acting  on  the  j-th  mode. 


respectively;  and  ftj  is  the  effective  frequency  of  the  j-th  mode 
related  to  the  fundamental  frequency  (ftj)  of  isolated  (A+B)  sub- 


system,  the  anharmonicity  (K*)  and  the  bath-induced  frequency  shift 
(6«j)  by 


IVJl-K'AWai  , 

*  a  J  o 


(  6) 


where  A^  is  the  steady-state  amplitude  of  the  j-th  mode  and  is  pro- 
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portional  to  the  applied  laser  field  intensity.  ' 

The  above  GLE  with  memory  effects  enables  us  to  study  the 
dynamics  of  LSSP  via  the  energy  relaxation  [governed  by  the 
damping  kernel  £3j(t)]  and  the  anharmonic  mode-mode  interaction 
function  N^j(t).  We  note  that  the  damping  kernel  characterizes 
the  dynamics  of  the  energy  relaxation  of  the  excited  mode  and  the 
energy  feedback  from  the  bath  modes,  while  the  dephasing  kernel 
characterizes  the  bath-induced  dephasing  of  the  excited  mode  in 
which  the  energy  population  does  not  change.  The  major  mechanisms 
which  cause  the  T^  and  T2  relaxations  are:  (i)  vibration-induced 
lattice-site  transitions  of  the  adspecies;  (ii)  anharmonic  mode¬ 
mode  coupling  of  the  (A+B)  subsystem  and  its  coupling  with  the 
bath  modes;  (iii)  the  interspecies  indirect  interactions (via  surface- 
phonon  modes) ;  (iv)  the  fluctuation  of  the  effective  dipole  of  the 
active  mode?  and  (v)  the  amplitude  and  direction  fluctuation  of  the 
effective  electric  field  acting  on  the  normal  coordinate  of  the 
system  caused  by  the  electron  charge  transfer  and  the  libration 
of  the  adspecies.  [For  a  more  detailed  discussion  on  the  overall 
line  broadening,  see  Ref.  3.3 


III .  Temperature-Dependent  Energy  Absorption  Rate 

Here  we  apply  the  GLE  described  in  the  previous  section  to 
the  calculation  of  the  energy  absorption  rate  of  the  system.  Two 
model  systems  are  considered:  (1)  there  is  only  one  active  mode  in 


the  (A+B)  subsystem  which  coupled  to  the  bath  modes,  and  (2) 
there  are  two  coupled  active  modes  in  the  (A+B)  subsystems  and 
both  are  coupled  to  the  bath  modes . 


By  linear  response  theory  and  time-dependent  perturbation 
theory,  we  may  express  the  temperature-dependent  averaged  energy 
absorption  rate  of  the  system  (via  the  active  dipoles)  as  follows 
(see  the  Appendix  for  the  derivation) : 

( <  -  e-'5H 

—i, , 


/  — \  = 


\  2too 


(7) 


where  F^j [m]  is  the  Laplace-Fourier  transform  of  the  velocity  cor¬ 
relation  function  defined  by 


F  .[»>:  =  ( 

rV  1 


**  -icot 


c 


(?A ' 


Thus,  we  see  that  for  a  system  with  more  than  one  active  mode, 
the  total  absorption  line  shape  is  a  linear  combination  of  the 
individual  line  shapes  (F^fu)])  plus  the  interference  term  (F^fw], 
i/  j),  characterized  by  the  Laplace-Fourier  transforms  of  the 
velocity  correlation  functions  F.  .(t)  which  will  be  computed  from 


the  GLE  [Eq. (5) ] . 

A.  Single  Active  Mode  System 

The  GLE  {Eq. (5) ]  in  this  single  active  mode  case  reduces  to 


Qfi )  *  fWj') Qtfxk'  +  j  M(i-i)  f «  >  +  f« )  ,  (9) 


%  i 


and  the  averaged  energy  absorption  rate  [Eq.(7)]  is  simply 


AFy  = 


We  now  solve  the  velocity  autocorrelation  function  F(t) 


(10) 

=  <Q(t)Q(0) > 


1 
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which  is  governed  by  the  GLE  in  the  absence  of  the  laser  field. 
By  multiplying  both  sides  of  Eq.(9)  by  Q(0)  and  performing  the 

ensemble  average  over  the  initial  conditions,  we  obtain 

•  ,t  ,t 

Fv'<  +  \  0.  00 
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Here  we  have  made  use  of  the  fluctuation-dissipation  theorem 
and  assumed  a  white-noise  process : 


F(0)=kT0/m,  < Q(0) fR(t) >  =0.  (12) 

By  taking  the  Laplace-Fourier  transform  of  Eq.(12),  we  find 


where  3  [u>]  and  M[w]  are  the  Laplace-Fourier  transforms  of  the 
damping  kernel  3(t)  and  the  dephasing  kernel  M(t) ,  respectively. 
In  obtaining  Eq. (13) ,  we  have  used  the  equilibrium  average, 

F (0)  =  <Q(0)Q(0)>  =  kTQ/m,  and  the  convolution  theorem  for  the 
Laplace-Fourier  transform. 

For  given  forms  of  the  memory  functions  /S(t)  and  M(t)  ,  we 
shall  be  able  to  find  the  related  Laplace-Fourier  transform  and 
in  turn  compute  the  averaged  energy  absorption  rate  by  the  real 
part  of  F[io],  We  note  that  in  using  the  convolution  theorem,  we 
require  the  memory  functions  6(t)  and  M(t)  and  the  velocity  auto¬ 
correlation  function  F(t)  to  be  well  behaved,  i.e.,  3(t),  M(t) , 

F ( t)  +  0,  as  t  +  “.  The  damping  kernel  which  governs  the 
energy  relaxation  and  affects  the  energy  population  of  the  active 
mode  is  chosen  in  the  form1"3 

/S(i)  =jS0  e  +  j-rcSU*  *(*)£)  J  ,  (l4~) 

where  3q  is  the  initial  damping  kernel  with  the  decay  constant  T 


■  ..  . 


j 

I 


'1 

J 


—2  2  1/2 

and  the  oscillation  frequency  =  [fi  -(T/2)  ]  '  .  Thus  damping 

kernel  satisfies  the  equation 

fid)  +  r/Sir  f  -  o  fr' 

for  the  underdamping  case  (i.e.,  T<  2 ST)  with  the  boundary  condi¬ 
tions  6(0)  =  6q  and  6 (0)  =  0 . 

We  note  that  this  oscillating-type  damping  kernel  contains 

the  memory  effects  of  the  T^  energy  relaxation  processes,  where 

the  excitation  and  relaxation  processes  of  the  active  mode  are 

accompanied  by  the  energy  feedback  of  the  bath  modes .  Such  memory 

effects  may  be  significantly  important  when  the  bath  modes  are 

only  slightly  excited  at  the  initial  stage  of  the  multiphonon 

relaxation  process  where  the  correlation  time  of  the  bath  modes 
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is  comparable  to  the  radiative  rate.  When  the  bath  modes  are 
highly  "heated"  due  to  the  rapid  energy  randomization  among  these 
modes  where  the  correlation  time  is  much  faster  than  all  the  other 
time  scales  of  the  processes,  the  damping  kernel  6(t)  loses  all 
memory  effects  and  the  energy  relaxation  processes  are  irrever¬ 
sible.  This  is  the  case  of  Markovian  processes,  whereby  the 
damping  kernel  (Eq.(14)]  reduces  to  a  Dirac  delta  function 

pi*:  =  t,'  Sot)  ,  <«•« 

(,  =  (fl/r) ,  (li.-i 

ij-’ 

where  is  a  constant  damping  factor  which  may  be  related  to  the 
microscopic  multiphonon  coupling  strength  and  the  Boson  occupation 
number  of  the  bath  modes  (see  Eq.(34)]. 

The  T 2  dephasing  kernel  M(t)  which  destroys  the  coherent 


nature  of  the  excitation  (i.e.,  the  stochastic  modulation  of  the 
fluctuating  frequency  of  the  active  mode) ,  but  does  not  change 

S  ,r 

the  energy  population,  is  chosen  in  an  exponential  form'' 

Mir  =  (-Q  t)  07 ' 

where  Y2  is  the  dephasing-induced  line  broadening  related  to  the 
phase  coherent  time,  tc,  by  Y2=Tc^’  Working  out  the  Laplace- 
Fourier  transforms  of  the  memory  functions,  A  [o> ]  and  M[u>],  and 
substituting  the  results  into  Eq.(13),  we  obtain,  from  the  real 
part  of  F[w],  the  temperature-dependent  averaged  energy  absorption 
rate  [Eq .  (10)  ]  . 

=(?E)V. 

where 


r(; v 1  &  h 

(i«) 

(-a) 

C  =  r<*>(V+XiU>)+ , 

M) 

X  =  27  +  {rAf , 

fac) 

1  = 

This  reflects  a  rather  general  line  shape  of  the  system,  where 
memory  effects  of  the  energy  relaxation  and  the  partially  coherent 
nature  of  the  excitation  processes  are  characterized  by  the  damping 
kernel  parameters  Bq  and  T  and  the  dephasing  rate  Y2* 

For  tractable  results  of  the  overall  line  broadening,  let 
us  consider  the  Markovian  limiting  case,  i.e.,  B  (t)  =  Bg<5(t).  The 
averaged  energy  absorption  rate  Eq.(18)  reduces  to  the  simple  form 


r  +  o)(i /+$p)S  " 


> 


\k  +  ft)  6  1 

A2  +  B*  J 


PO 


e>  =  ^cr;-r2).  f«.<0 

The  important  features  of  the  above  temperature-dependent  line 

shape  are:  (1)  for  low  temperature  (kTQ<<ftcj),  P(tq)  is  independent 

to  Tq  when  the  overall  broadening  ^  +  y  )  is  weakly  temperature- 

dependent  [in  general,  is  temperature-dependent  -  see  Eq.(34)]; 

(2)  the  asymmetric  line  shape  is  due  to  the  dephasing  term  y  A 
which  is  antisymmetric  with  respect  to  ai=Q+ (y^y2)  /  (2f2)  ,  i.e., 

A  >  0  when  u><u,  and  A<0  when  a>>tu ?  and  (3)  for 

YlY2<<a>'  Ec3  -(21)  reduces  to  essentially  a  Lorentzian  with  the 
FWHM  =  Yi+Y2»  where  for  low  excitations  (coherent  processes),  the 
overall  line  broadening  is  dominant  by  y^,  while  for  high  excita¬ 
tions  (incoherent  processes)  the  strong  interactions  between  the 
active  mode  and  the  dense  bath  modes  destroy  the  coherent  nature 
of  the  excitations,  and  the  T2  dephasing  broadening  (y2)  is  domi¬ 
nant.  We  note  that  the  simple  additive  nature  of  the  overall 
broadening  [the  (y1+y2)  term  in  Eq.(2f)]  is  available  only  for 
the  Markovian  processes,  where  the  correlation  time  of  the  bath 
modes  is  much  faster  than  the  and  T2  relaxation  rates  and  the 
memory  effects  of  the  damping  kernel  are  not  significant.  We 
believe  that  the  dynamical  transition  of  the  nature  of  LSSP  (from 
a  coherent  to  an  incoherent  excitation)  may  be  described  by  the 
relative  magnitudes  of  the  relaxation  parameters  Bg,  T  and  y2 
contained  in  the  line  shape  [Eqs.(18)  and  (21 ) 1 >  which  reflects 
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not  only  the  internal  level  structure  of  the  active  mode  (which 
is  characterized  by  the  conformation  of  the  adspecies  and  the 
heterogeneity  of  adsorbent)  but  also  the  dynamical  behavior  of 
the  excited  adspecies,  e.g.,  vibration-induced  surface  migration, 
surface  recombination  and  the  rate  processes.^ 

B .  Two-Coupled  Active  Modes 

For  the  situation  that  there  are  two  active  modes  with  fre¬ 
quencies  close  together  (within  the  detuning  range  of  the  field) , 
then  both  modes  will  strongly  interact  with  the  field.  While 
acting  with  each  other  (via  linear  coupling) ,  these  two  modes  both 
are  coupled  to  the  bath  modes  (via  multiphonon  processes) ,  and  the 
energy  deposited  in  the  active  modes  will  gradually  relax  to  the 
bath  modes.  At  steady-state,  the  energy  stored  in  the  active 
modes  remains  unchanged  and  the  averaged  energy  absorption  rate  of 
the  system  must  equal  the  energy  dissipative  rate  of  the  active 
modes.  We  may  express  this  steady-state  averaged  power  absorption 
classically  by"^ 


where  f^(t)  is  the  generalized  force  for  the  normal  coordinate 
governed  by  the  usual  Langevin  equation  (without  including  the 
memory  effects) ;  p^  and  P2  are  the  averaged  power  absorption  for 
the  uncoupled  active  modes  and  p^  is  the  interference  term  due 
to  their  couplings. 

For  a  quantum  system  with  memory  effects,  instead  of  Eq.(23), 
the  (time  and  ensemble)  averaged  power  absorption  is  given  by 
[from  Eq .  (7) ] 
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<jj>  =  <l>  -  <t>  *  <%? +  <t> ,  (m 

<%>  =  W  ( %  \ E; Ej) ( fe  RjM)  ,  ^ 

which  is  similar  to  that  of  the  classical  expression,  except  that 
here  we  take  the  ensemble  average  <•••>  over  the  quantum  states, 
and  the  instantaneous  velocity  Q(t)  is  replaced  by  its  correlation 
function  ^ (t)  = <Qi (t) (0) >  due  to  the  memory  kernels. 

To  obtain  the  above  quantity,  we  first  consider  the  GLE 
CEq - ( 5 ) ]  for  the  case  of  two  coupled  active  modes: 

Q(i  *  +  S2*\  \  =  +  //n,  , 

°  »  4 

iajtu  fya-t,3/> A'-rlj 

6  6  i  ^  # 

By  multiplying  both  sides  of  the  above  equations  by  (^(0)  and  Q2(0) 

and  performing  the  ensemble  average,  we  obtain  a  system  of  four 

equations  for  the  velocity  correlation  functions  F^tt),  F22(t), 

F12(t)  and  F21(t),  which  in  turn,  by  taking  the  Laplace-Fourier 


transform,  give  us  the  coupled  equations 


-r  j  [«•  .1  +  S2*  N )  F,  ~  ~  f  /  t 

(  ^  *-  +  -T£ Mia 0  ~  4  f\uCuj;  p,,iwi  -  i  T0  h?x  ,  UU) 

\  '  r  F,^’  ■  Kal«i  F21‘:ti  =  0,  (>:•<  > 

(  i v  r  i  Fj^l  “  ^  A,j.[^l  fv~i  =  o. 

Again,  we  have  used  the  fluctuation  theorem  and  assumed  a  white- 
noise  :  F^j  (0)=<r  jkT0/rni  and  <  Q^O)  f  J(t)>=0. 


(JJ.C) 


('■:-<  > 


The  solutions  of  Eq.(28)  are  straightforward,  where  by 
Cramer's  rule  we  obtain 


C,/('n,V)7 

(?7'0 

F>:  -  !>»:  =  4%  l(x»:/p  , 

fif.C  > 

V  =  //•’,  n  c1;  -  NS*:  , 

) 

C!,  =  -  eru>~  +  • 

t'f.e / 

In  principle,  for  given  forms  of  the  memory  functions  8^(t), 

M^(t)  and  N^2(t)  we  shall  be  able  to  compute  their  Laplace-Fourier 
transforms  and  thus  obtain  the  line  shape  function. 

For  tractable  results,  let  us  investigate  the  Markovian  pro¬ 
cesses  and  the  memory  functions  with  the  following  forms: 


j}Aii )-&;£(*>  * 

(  jJ. 

M;  >  •  =  to'f'C-Xi*)  , 

A, a  } 

where  8^  and  y^  represent  the  line  broadening  of  the  active  modes 
achieved  by  the  T1  (energy)  and  T2  (phase)  relaxation,  respec¬ 
tively,  and  y12  is  the  reciprocal  correlation  time  of  the  inter¬ 
action  between  these  two  active  modes  with  the  mode-mode  coupling 
strength  2  related  to  the  interaction  potential  V(Q^,Q2)  by 
2X^2 =  O  V/3Q19Q2)q.  Substituting  the  Laplace-Fourier  transforms 
of  the  above  memory  functions  into  Eq.(29),  the  real  parts  of  the 
results  give  us  the  averaged  power  absorption  [from  Eq . (24) ] 
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<j|>  =  P(T.l  l  If(“)  , 

re  OC^r^KO2.  +  +  i&tKOz* 

Ij2m=  -  > 

k 

Cr^j  ~  RA3j~  12  ^i  j 

Ka-  +  Ay  , 

2^2=  f<L(2  +  2°^XlL'2,\  f 

L,  =  AAz'&A  , 

L2  -  ^,5*  +  /l2  £■)  ,  2 

K,^(i,MrcS A^A*  , 
13,^-  as  f  +  /s,.»  ( r,s  \>  V  ^  +  A.>  /».,.  ] , 

R  -  t*- . 
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(32) 
(33. «) 

Gitf) 

(33.  c) 

(ii.e) 
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(»i) 


The  above  line  shape,  while  complicated,  is  a  rather  general 
response  function  of  a  system  with  two  coupled  active  modes  and 
enables  us  to  describe  the  dynamical  nature  of  LSSP,  where  both 
the  memory  effects  of  the  dephasing  kernel  (governed  by  Yj  and  y2) 
and  the  mode-mode  interaction  (governed  by  Y]^  are  included. 

Using  these  expressions  [Eqs . ( 31) - (33) ] ,  we  may  generate  different 
line  shapes  by  varying  the  relative  magnitudes  of  the  parameters 
y^»  Y^2  an<*  *12*  s*nce  the  experimental  results  of  LSSP  are 
not  available  to  date,  we  shall  not  show  here  the  numerical  plots 
of  the  line  shapes  for  an  actual  system.  However,  some  important 
features  are  pointed  out  as  follows:  (i)  the  line  shape  described 
by  Eq.(31)  is  essentially  composed  of  two  Lorentzian-type  line 


shapes  (with  the  asymmetry  due  to  the  dephasing  where 

the  peak  values  are  located  near  the  frequencies  of  the  active 
modes;  (ii)  the  overall  line  shape  is  not  only  characterized  by 
the  coupling  strength  X^  (which  governs  the  interference  term  in 
but  also  by  the  frequency  separation  of  these  two  active 
modes  -  51^);  and  (iii)  for  the  very  weak  coupling  case  (X^'  °» 
due  to  the  big  separation  between  the  frequencies  of  the  active 
modes)  with  the  laser  frequency  near  resonant  to  one  of  the  active 
mode  frequency,  Eq. (31)  reduces  to  the  single-mode  line  shape 
[Eq .(21)1. 

IV.  Multiphonon  Processes  and  Laser-Selective  Effects 

In  this  section  we  shall  present  a  classical  model  to  study 
the  nature  of  the  laser-selective  effects  via  the  multiphonon 
processes  which  have  been  recently  developed  by  quantum  models. 

For  an  N-atom  polyatomic  molecule  chemisorbed  on  a  solid  surface, 
the  number  of  "frustrated"  vibrational  normal  modes  is  3N.  We  now 
divide  these  3N  intramolecular  modes  into  two  types:  the  active 
(A)  mode(s)  and  the  inactive  surface-phonon  (B)  modes.  The  bulk- 
phonon  modes  with  lower  frequencies  are  considered  to  provide  the 
bath  (C)  modes  for  the  (A+B)  modes.  In  this  picture,  the  laser 
photon  energy  is  first  deposited  into  the  active  (A)  mode,  and 
gradually  leaks  from  this  mode  to  the  B  and  C  modes  accompanied 
by  some  energy  feedback  from  the  B  modes.  We  then  may  define  the 
laser-selective  excitation  (or  local  heating)  in  a  nonequilibrium 
fashion,  where  the  selectively  excited  A  mode  (or  the  local  heating 
of  the  B  modes)  store  the  photon  energy  in  the  specific  vibrational 
degrees  of  freedom  on  a  time  scale  long  enough  for  the  related 


t 


chemical  reaction  or  the  bond  breaking  to  occur. 

56 

In  previous  papers,  '  we  have  discussed  this  excitation 

nature  (selective,  local  heating  and  nonselective  thermal  heating) 

by  the  dynamical  energy  populations  of  the  system.  It  was  shown 

that  the  laser-selective  effects  are  characterized  by  the  pumping 

rate  and  the  energy  relaxation  rates  of  the  A  and  B  modes,  where 

the  intramolecular  vibrational  relaxation  (IVR)  plays  an  important 

role  in  the  selective  excitation  of  the  A  mode  and  local  heating 

of  the  B  modes.  The  vibrational  level  width  induced  by  IVR  via 

8  17  18 

multiphonon  processes  is  given  by  '  ' 

(}<u> 

*  i 

irw>  •  (m 

i  a  6 

<n^>  is  the  thermally-averaged  Boson  occupation  number  of  the  j-th 
mode  (with  frequency  )  which  is  strongly  coupled  to  the  active 
mode;  measures  the  strength  of  the  multiphonon  coupling  and  is 
proportional  to  the  (p+l)-th  derivative  of  the  intramolecular  an- 
harmonic  coupling  potential  (for  a  p-order  multiphonon  process) 
evaluated  at  the  equilibrium  point;  and  the  Dirac  delta  function 
simply  implies  the  conservation  of  energy  as  in  the  Fermi  golden 
rule . 

We  now  find,  by  a  fully  classical  model,  the  required  condi¬ 
tion  for  a  slow  IVR  rate  (i.e.,  a  small  damping  y) ,  which  in  turn 
gives  us  the  possibility  of  selective  bond  breaking.  We  begin  with 
the  nonlinear  coupled  equations  of  motion  of  the  normal  coordinates 
Qa  (for  the  active  mode)  and  Qj  (for  the  bath  modes): 
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(ir.4) 


X(p)  is  the  p-phonon  coupling  constant  related  to  the  derivative 
of  the  mode-mode  interaction  potential,  V,  by 

f=bM{^'vAdQAT3Qh]D ,  <**) 

r  ° 

and  p  =  It  (&  .=integer)  is  the  order  of  the  multiphonon  coupling, 

j  3  3 

i.e.,  the  active-mode  energy  gradually  leaks  into  the  bath  modes 

by  exciting  p  phonons  phonons  of  the  j-th  mode  with  frequency 

fi.,  etc.).  With  the  boundary  conditions  Q,  (0)  Q*  (0)  ^  0 ,  Q.(0)=A. 

3  A  A  3  3 

and  6^(0)  =  0,  we  shall  find  the  condition  (or  the  relation  among 
the  frequencies  of  the  intramolecular  modes)  for  the  energy  initially 
stored  in  the  active  (A)  mode  to  be  significantly  transferred  to  the 
bath  modes.  The  coupled  equations  may  be  solved  by  an  iterative 
scheme  provided  the  multiphonon  coupling  may  be  treated  as  a  small 
perturbation.  Note  that  the  perturbative  method  is  usually  valid 
if  one  considers  a  high-order  multiphonon  process  where  the  coupl¬ 
ing  strength  X ^  is  a  strongly  decreasing  function  of  p.7'8  For 
the  strong  coupling  case,  e.g.,  (linear)  single-phonon  process,  we 
must  solve  the  coupled  equations  exactly. 

By  substituting  the  zero-th  order  solution  (t)  =  A. cos  (fi  .  t) 

J  J  J 

into  Eq.(35.a),  we  obtain  the  first-order  solution  for  QA(t) : 


17 


q%  =  tyo)  cm. (n^i)  +  [ ^ +  Td) ,  Mo 


iin  =  Am  yt'  J Af  . 


o  a 


ru.tH) 


The  integral  I(t)«>cX'p  ,  characterizing  the  strength  of  the  multi¬ 
phonon  coupling,  is  quite  complicated  for  arbitrary  values  of  the 
integers  .  For  simplicity  and  without  losing  the  spirit  of  our 
resonant  condition,  let  us  consider  the  p-phonon  process  with 
^1  ~  ^2  *  *  * an<^  ^l=  *"2  *  *=  ^  j  *  T^e  couPHn9  integral  then  appears  as 

f  x,f) 

}  — ; — 7 - T^r  /  (37) 

which  clearly  indicates  that  as  X  p  -*-0 ,  then  I(t)-*-0,  unless  the 

exact  "internal  resonant"  condition  is  fulfilled,  i.e.,  ft.  =  |p-2k|ft_. 

a  J 

(where  p  and  k  are  integers) .  However,  when  the  coupling  strength 

X(p)  has  a  small  finite  value,  the  coupling  integral  I(t)  will  be 

significant  when  ft.  is  near  "internal  resonant"  to  ft.,  i.e., 

A  J 

ftA= | p— 2k | ft  j .  Thus  the  necessary  condition  of  the  IVR  to  be  signifi¬ 
cant  is  the  near  "internal  resonant"  condition,  which  is  somehow  a 
weaker  condition  than,  that  of  the  exact  resonant  condition  governed 
by  the  delta  function  [ Eq . (34) ] .  In  fact,  if  we  approximate  the 
upper  time  of  the  coupling  integral  [Eq.(36.a)l  from  zero  to  infinite,  we 
obtain  the  delta  function  6 (ft. - (p-2k) ft  .  ] . 

A  j 

Therefore,  the  rate  of  IVR  is  restricted  by  two  factors: 

(i)  a  low  order  of  the  multi-quantum  processes  and  (ii)  the  near 
"internal  resonant"  condition,  i.e.,  ftA-E)ljftj  s0.  This  means  that 
a  rapid  IVR  rate  (i.e.,  a  large  line  broadening  which  reduces  the 
energy  absorption  rate)  may  be  achieved  only  when  both  of  the  above 


Pi** mf  - 


conditions  are  fulfilled;  otherwise,  we  expect  a  much  slower  IVR 
rate  which  in  turn  cuts  off  the  line  broadening  and  makes  the 
selective  bond  breaking  possible. 

V  Discussion  and  Conclusion 

In  our  previous  work,^'^  we  have  investigated  the  selective 
versus  nonselective  nature  of  LSSP  in  terms  of  parameters  such 
as  the  pumping  rate,  coupling  factor  and  damping  rate.  We  shall 
now  present  some  possible  physical  systems  exhibiting  selective 
effects  and  surface-enhanced  bond  breaking. 

Let  us  consider  a  system  which  may  be  pictured  as  a  set  of 
groups  in  the  frequency  domain.  For  simplicity,  we  shall  consider 
the  case  of  only  two  groups,  A  and  B,  in  which  group  A  is  referred 
to  as  the  excited  broup  which  consists  of  the  active  mode  plus 
some  intramodes  of  the  adspecies  coupled  strongly  to  the  pumped 
mode,  and  group  B  consists  of  the  remaining  modes  of  the  total 
system  (adspecies  plus  solid) .  Depending  on  the  relative  magni¬ 
tudes  of  the  intragroup  (Rlntra)  and  the  inter  group  (Rlnter) 
coupling  rates,  which  characterize  the  energy  randomization  rate 
wi thin  and  between  the  group(s) ,  respectively,  and  the  laser  pump¬ 
ing  rate  (V),  we  may  introduce  several  types  of  laser  excitations: 
(i)  mode  (bond) -selective ,  (ii)  molecule-selective  and  (iii)  non¬ 
selective  thermal  processes. 

For  LSSP  to  be  characterized  as  a  selective  type  (mode-  or 
molecule-selective) ,  we  shall  require  the  pumping  rate  to  be  faster 
than  the  energy  relaxation  rate(s),  i.e.,  V >  R*ntra  and/or 
V  >Rinter  >>Rintra.  For  the  case  of  very  high  (low)  pumping  (re¬ 
laxation)  rate,  V»Rintra,  Rinter,  a  mode-selective  type  excitation 
would  be  possible.  However,  for  low  intensity  excitations,  the 


mode-selective  seems  to  be  less  likely  than  the  molecule-selective 
where  only  the  condition  required  is  V,  Rin^ra>>  Rln^er.  By  the 
concept  of  the  "energy-gap  law,"  it  is  highly  possible  for  a  system 
to  have  very  weak  intergroup  coupling  but  strong  intragroup  coupl¬ 
ing,  e.g.,  a  large  molecule  adsorbed  on  a  solid  surface  like 
SF^-metal  and  a  long  chain  adspecies  like  A-B-C-D-metal .  In  these 
systems  the  vibrational  stretch  of  the  excited  species  (which  is 
not  directly  connected  to  the  surface)  may  behave  like  a  small 
subsystem,  where  the  excited  bond  coordinate  is  composed  by  the 
normal  coordinates  of  the  related  group.  This  selectively-excited 
adspecies,  while  strongly  coupled  within  the  group,  is  weakly 
coupled  to  the  remaining  species  (or  functional  group)  of  the 
system.  Therefore,  selective  bond  breaking  of  adspecies  (or  cer¬ 
tain  functional  group)  will  be  possible. 

Another  important  feature  of  selective  bond  breaking  of  a 
heterogeneous  system  (and  usually  not  present  in  a  gas-phase  homo¬ 
geneous)  is  the  surface-enhanced  local  field  acting  on  the  adspecies. 
The  pumping  rate  could  be  actually  enhanced  by  a  significant  factor 
when  the  local  field  is  increased  due  to  surface  effects  such  as 
roughness  and  electron  transfer  of  the  substrate.  By  the  simple 

relation  of  the  pumping  rate  and  the  local  electric  field 
2 

I  *  (E^qc)  •  we  are  at>le  to  see  that  the  required  laser  intensity 

Q 

for  selective  excitation  may  be  reduced  to  a  factor  of  10  when 

4 

the  local  electric  field  is  enhanced  by  a  factor  of  10  .  We  note 
that  the  enhancement  of  the  local  surface  field  which  governs  the 
experimental  evidences  of  surface-enhanced  Raman  scattering  shall 


also  play  an  essential  role  in  LSSP,  where  considerably  lower  laser 

3  2 

intensities  (10-10  W/cm  )  are  used  as  compared  with  that  of  photo¬ 
dissociation  of  polyatomic  molecule. 

In  conclusion,  we  should  mention  that  selective  laser-stimulated 
desorption  (bond  breaking)  processes  are  characterized  not  only  by 
the  coherent  nature  of  the  laser  radiation  but  also  by  the  hetero¬ 
geneity  of  the  substrate  (e.g.,  roughness,  orientation  and  charge 
transfer)  and,  very  likely,  by  the  internal  excitation  structure  of 
the  pumped  state  (e.g.,  the  miqration-induced  and  vibration-trans¬ 
lation-coupled  quasi-continuum  state  near  the  desorption  region) . 

We  finally  note  that  the  excitation  nature  of  LSSP  is 
characterized  by  the  dominant  one  (or  more)  component  of  the 
above  (i)  -  fc'fi)  processes,  which  in  turn  is  governed  by  the 
relative  magnitudes  of  the  pumping  rate,  the  and  T2  relaxa¬ 
tion  rate  and  the  dynamical  behavior  of  the  excited  modes  (via 
the  memory  effects  of  GLE) . 


Consider  a  system  of  N  IR  active  modes  with  the  initial 
quantum  state  |k>,  which  interacts  with  an  electric  field  of 
frequency  w,  resulting  in  transitions  to  some  final  quantum 
state  |k’>.  Let  the  laser  radiation  be  linearly  polarized  in  the 
z  direction,  whereby  we  write  the  laser  monochromatic  field  as 

=  /.l) 

where  Eq  is  the  amplitude  of  the  field.  We  assume  here  that  the 
field  is  spatially  uniform,  since  the  wavelength  is  large  com¬ 
pared  to  the  dimensions  of  the  adspecies.  The  interaction  Hamil¬ 
tonian  between  the  field  and  the  adspecies  can  be  written  as 

H'lto  -  £  Mi  E-  }  (A.2) 

A 

where  E^  is  the  effective  amplitude  of  the  electric  field  acting 
on  the  i-th  normal  modes  with  the  active  dipole  moment  operator  vk  . 

According  to  time-dependent  perturbation  theory,  the  transi¬ 
tion  probability  of  a  transition  from  the  state  |k>  to  the  state 
1 9 

| k  *  >  is  given  by 

where  i-^'  an<^  the  laser  field  is  turned  on  slowly  with 

an  exp(nt)  factor  (with  n-*-0+)  .  Using  the  expression  for  the  inter¬ 
action  Hamiltonian  H'(t)  [Eq.(A.2)3  and  working  out  the  integral 

of  Eq . (A. 3) ,  we  obtain  the  transition  rate  (i.e.,  the  transition 
probability  per  unit  time)  as 
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which  is  still  time  dependent.  However,  we  shall  be  interested  in 
the  time-averaged  quantity  (over  the  field  cycle  T  =  2tt/w)  .  Using 
the  fact  that  the  time  average  of  the  terms  sin(2ut)  and  cos  (wt) 
vanish  and  transforming  the  rhs  of  Eq.(A.4)  into  the  Dirac  delta 
function , 

[;  m  J (-IL* u)  ,  M-s) 

we  obtain  the  time-averaged  transition  rate  (which  is  now  time 


independent)  as 
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We  then  may  calculate  the  time-averaged  energy  absorption 
rate  of  the  system  given  by 


<77>  = 
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where  p^  is  the  probability  that  the  system  was  in  the  initial  state 
|k>.  Since  the  summations  over  k  and  k’  range  over  all  the  quantum 
states  of  the  system,  we  may  interchange  these  indices  in  the  sum¬ 
mation  over  the  second  delta  function  in  Eq.(A.5)  to  obtain 


a  < 


By  assuming  the  system  initially  to  be  in  thermal  equilibrium,  we 
can  relate  the  thermal  probabilities  of  the  system  p^  and  p^ ,  by 
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tc  =  ,  & V 

where  B  =  (kTfl )  ^  and  Tg  is  the  initial  temperature  of  the  system. 
In  order  to  convert  Eq.(A.8)  to  the  Heisenberg  picture,  we  intro¬ 
duce  the  Fourier  transform  of  the  delta  function, 

_c>o 

Using  the  expressions  of  Eqs.(A.9)  and  (A. 10) ,  Eq.(A,8)  becomes 
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Since  the  |k>  and  | k  * >  are  eigenstates  of  the  unperturbed  system 
[i.e. ,  H' (t)  =  0] ,  we  have 

-xOJ 


"  |^>=  [>)  K>  , 

(/!  =  0'\  UW  > 


(AAVi 


where  U  (t)  =  exp (-iH0t/ft)  is  a  unitary  operator  of  the  unperturbed 
Hamiltonian  HQ .  Combining  Eqs.(A.ll),  (A. 12)  and  (A. 13)  and  using 
closure  relation  of  the  final  state  j  k  * > ,  we  obtain 

{^f>  -  (l  - e^u)  \  H  i:mt<  fyitmi')  .  M.W') 

l4  0 

where  vu  (t)  is  the  dipole  moment  of  the  i-th  active  mode  defined 
in  the  Heisenberg  picture  by 

Mid)  =  Life  M/o){J(t)  (A-'r.  / 

A 


<•••>  denotes  the  ensemble  average,  i.e.,  in  writing  Eq.(A.14), 
we  have  used  the  equilibrium  ensemble  average 


}'  %  (4lAtt)M(0)lt>  =  CM(ir)Mwy  .  (AJ6) 


Thus,  we  have  shown  that  the  time-dependent  (time  and  ensemble) 
averaged  energy  absorption  rate  [Eq.(A.14)l  is  characterized  by 
the  initial  temperature  Tq  of  the  system  and  the  Fourier  transforms 
of  the  time-correlation  function  of  the  dipole-moments  of  the 
absorbing  adspecies  in  the  absence  of  the  laser  field. 

In  order  to  employ  GLE  which  gives  the  correlation  function 
of  the  velocity,  we  shall  further  express  the  averaged  absorption 
rate  in  terms  of  the  velocity  correlation.  To  do  this,  let  us 
expand  the  dipole  moment  in  terms  of  the  related  normal  coordinates 


^  (*)  + 


which  gives  us 


gf. <<?•«>$*'>, 

0  *  J 


i a.i ?; 


where  q^ ,  q^  are  the  effective  classical  charges.  Here  we  have 

used  q^  =  ,  q j  .  =  Uj0  and  dropped  the  yio  ,  terms 

which  are  related  to  the  pure  rotational  transitions  and  give  no 

contributions  to  the  vibrational  excitation  processes  (by  the 

selection  rule).  Substituting  of  Eq. (A. 18)  into  Eq.(A.14)  and 

using  the  relation  between  the  displacement  correlation  and  the 

.  .  2 

velocity  correlation  <Q(t)Q(0)>  =  u  <Q(t)Q(0)>,  we  obtain 

where  R*F^j[u)  is  the  real  part  of  the  Laplace-Fourier  transform  of 
the  velocity  correlation  function  defined  by 

Fxjt“l=  pt  <4(«q,(.)>.  (A.io'. 

We  note  that  the  velocity  correlation  function  is  even  in  time, 
i.e.,  FAj(-t)  -  FAj(t) . 
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<•••>  denotes  the  ensemble  average,  i.e., 


in  writing  Eg. (A. 14) 


we  have  used  the  equilibrium  ensemble  average 

I  %  <4lMltW0)U>  =  <Mlt)M'0)y  . 


(A.) 6) 


the  velocity  correlation  function  defined  by 

t)  ,  fyi)  =  <  Q/t)  Q3<o)>  .  (A. 20] 

We  note  that  the  velocity  correlation  function  is  even  in  time 
i.e. ,  F.  . (-t)  *  F. . (t) . 


